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Atomic nuclei and nano-scale metallic grains are in the crossover regime 
of pairing correlations between the bulk limit, where the Bardeen- 
Cooper-Schrieffer (BCS) theory of superconductivity is valid, and the 
fluctuation-dominated regime, where BCS theory breaks down. In 
this fluctuation-dominated regime, the pairing gap is comparable to or 
smaller than the single-particle mean level spacing. We discuss ther- 
mal signatures of pairing correlations in nuclei and ultra-small metallic 
grains that survive despite the large fluctuations of the pairing field. 



1. Introduction 

Pairing correlations lead to superconductivity in bulk metals. Effects of the 
pairing correlations in nuclei, such as a gap in the excitation spectrum of 
even-even nuclei, are well documented. Superconductivity was explained by 
the Bardeen-Cooper-Schrieffer (BCS) theoryP Following its introduction in 
electronic systems, the BCS approximation was applied to nuclei by Bohr, 
Mottelson and Pineal and by BelayevP 

Single-electron tunneling spectroscopy experiments in ultra-small metal- 
lic grains connected to external leads probed the discrete spectra of these 
grainsP A gap was identified in the excitation spectra of larger grains 
with an even number of electrons. This led to extensive studies of pair- 
ing correlations in nano-scale metallic grainsP Recent technical advances 
are providing better experimental control over the size and shape of these 
grains P 

BCS theory is a mean-field theory valid in the limit where the pairing 
gap A is much larger than the single-particle mean level spacing 5. How- 
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ever, in a finite-size system, fluctuations of the order parameter around 
its mean-field solution can be important. As the linear size of a metallic 
grain decreases, its single-particle mean level spacing S increases. Anderson 
arguecP that the smallest size at which a metallic grain remains a super- 
conductor corresponds to S ~ A. In the smallest metallic grains studied in 
Ref. HI the excitation spectrum of the even particle-number grain did not 
exhibit a noticeable gap when compared to the excitation spectrum of the 
grain with an odd particle number. These grains belong to the fluctuation- 
dominated regime A < 5, in which fluctuations of the pairing field become 
important and BCS theory breaks down. It was proposed that signatures 
of pairing correlations in this regime can still be identified through the 
particle-number parity dependence of thermodynamic observablesP^ 

In a nucleus, the gap is typically of the order of the single-particle mean 
level spacing or somewhat larger. We therefore expect fluctuations in the 
pairing field to also be important in nuclei. 

Here we discuss thermal signatures of pairing correlations in both nu- 
clei and ultra-small metallic grains. We use methods that go beyond the 
mean-field or BCS approximations and take into account both thermal and 
quantal fluctuations. These methods are briefly described in Sec. [2j Appli- 
cations to nuclei are presented in Sec. [3] In Sec. 13.1 I we discuss the effective 
configuration-interaction shell model Hamiltonian used in our studies, and 
in Sec. 13. 21 we review the auxiliary-field Monte Carlo (AFMC) method used 
to calculate thermal and statistical nuclear properties in very large model 
spaces. In Sees. 13. 3\ 13.41 and 13.51 we discuss, respectively, signatures of 
pairing correlations in the heat capacity, the spin distribution of nuclear 
levels and the thermal moment of inertia. In Sec. |4] we present applications 
to nano-scale metallic grains. In Sec. 14.11 we discuss the universal Hamilto- 
nian describing the effective low-energy Hamiltonian of grains whose single- 
electron dynamics is chaotic. Applications of AFMC to calculate signatures 
of pairing correlations in thermodynamical observables of metallic grains 
are discussed in Sec. 14.21 The single-particle Hamiltonian of a chaotic grain 
follows random-matrix theory, and physical observables of the grain un- 
dergo mesoscopic fluctuations. The study of these mesoscopic fluctuations 
requires calculations for a large number of realizations of the single-particle 
spectrum, and an efficient finite-temperature method is therefore needed. 
We discuss such a method in Sec. 14.31 combining spin and number-parity 
projections with thermal and small-amplitude quantal fluctuations of the 
pairing field. In Sec. l4.4l we use this method to study the mesoscopic fluctu- 
ations of the heat capacity and spin susceptibility of an ultra-small metallic 
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grain. We conclude in Sec. [5] by comparing thermal signatures of pairing 
correlations in nuclei with similar signatures in nano-scale metallic grains. 

2. Finite-temperature methods: beyond the mean field 

Correlations beyond the mean-field approximation at finite temperature can 
be taken into account systematically using the Hubbard- Stratonovich (HS) 
transformation,^ i n which interaction effects are described by including 
fluctuations around the mean-field solution. 

Hubbard- Stratonovich transformation. The equilibrium density matrix 
e —PH describing a system with an Hamiltonian H at temperature T = 1//3 
can be interpreted as the imaginary-time propagator with /3 playing the 
role of imaginary time. The HS transformation 



expresses this imaginary-time propagator as a coherent superposition of 
one-body propagators U a with a Gaussian weight G a . Each U a describes 
the imaginary-time propagator of non-interacting particles moving in time- 
dependent external auxiliary fields cr(r). We note that the HS decomposi- 
tion can be done either in the particle-hole channel or in the particle-particle 
channel. 

The HS decomposition (JTJ can be used to calculate various thermal 
observables. The partition function Z(T) = Tre~ H / T is found by taking 
its trace 



The mean-field approximation is obtained by evaluating the integral in 
([2]) in the saddle-point approximation. In the particle-particle decomposi- 
tion, this leads to the BCS approximation. 

Static path approximation. To go beyond a mean-field approximation, it 
is necessary to include fluctuations of the auxiliary fields a. At high temper- 
atures, it is sufficient to include static (thermal) fluctuations of the a fields, 
an approximation known as the static path approximation (SPA). EES Of 
particular importance are large-amplitude fluctuations of the order param- 
eter. For example, in the Landau theory of the nuclear shape transitions, 
static fluctuations in the quadrupole shape parameters were found to be im- 
portant for understanding the observed temperature and spin dependence 
of the giant dipole resonance.^ 
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Static path approximation plus random phase approximation. The SPA 
can be improved by including small-amplitude time-dependent (quantal) 
fluctuations around each static configuration crP^JEl] This can be accom- 
plished by expanding ct(t) = J2 r Vre lUrT where uj r = 2irr/(3 (r integer) are 
the bosonic Matsubara frequencies. For any given static value o~o, the in- 
tegration over ay (r 7^ 0) is carried out in the saddle-point approximation, 
resulting in an <7o-dependent random phase approximation (RPA) correc- 
tion factor. Finally the integration over the static o~q is carried out exactly 
(i.e., including large- amplitude static fluctuations). This approximation 
is known as the SPA+RPA. The most important fluctuations originate in 
those of the order parameters. The approximation breaks down below a cer- 
tain critical temperature under which the Gaussian fluctuation in a given 
o r (for r^O) becomes unstable. 

Auxiliary-field Monte Carlo (AFMC). To account for correlation effects 
in full, it is necessary to include all fluctuations - both thermal and quantal 
- of all the auxiliary fields a (including large-amplitude quantal fluctua- 
tions). This requires an integration over a very large number of a fields (at 
all time slices), and in practice can only be done by Monte Carlo methods. 
Such a quantum Monte Carlo method is generally known as the auxiliary- 
field Monte Carlo (AFMC) method and has been used in strongly corre- 
lated electron systems!^ In the context of the configuration-interaction 
shell model the method is known as the shell model Monte Carlo (SMMC) 
methodPHM] 

3. Nuclei 

3.1. Nuclear Hamiltonian 

Here we use the framework of the configuration-interaction shell model 
approach. The single-particle energies are derived from a central Woods- 
Saxon potential plus spin-orbit termP^l Our interaction includes dom- 
inant components^ of effective nuclear interactions: monopole pairing 
plus multipole-multipole interactions (quadrupole, octupole and hexade- 
cupole). 27 The latter are obtained by expanding the separable surface- 
peaked interaction u(r,r') = — x{dV / dr)(dV / dr')5(i — ?') (V is the central 
Woods-Saxon potential) into mulipoles. The coupling constant \ is deter- 
mined self-consistently^S 
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where p is the nuclear density. The quadrupole, octupole and hexadecupole 
interaction terms are retained and renormalized by factors of 2, 1.5 and 1, 
respectively. 

3.2. AFMC 

For finite-size systems it is often necessary to calculate observables at fixed 
particle number. The thermal expectation value of an observable O at fixed 
particle number A is given by 

= Tr A (Oe-f> H ) JP[#A(Q). 

where Tr^ denotes a trace at fixed particle number A and we have used the 
HS transformation (Q}. Here W a — G a \TTAU a \ is a positive-definite weight 
function, <& a = Tr^ U a /\TrAU a \ is the Monte Carlo "sign" and (0) c = 
TiA(OU a )/TTAUa- The sample-specific quantities TtaUg and {G) a can 
be evaluated using matrix algebra in the single-particle space. We denote 
by U CT the N sp x N sp matrix representing U a in the single-particle space 
containing N sp single-particle orbitals. In the grand-canonical ensemble, 
we then have 

TrUa- = det(l+U CT ) , (5) 

and 



(6) 



For a one-body operator O = J^ij {AO\j) a \ a ii the grand-canonical expec- 
tation value (0) a can be calculated using ((6]). For a two-body operator, 
Wick's theorem can be used to express the two-body expectation values in 
terms of one-body expectation values. 

Canonical quantities can be evaluated using particle-number projection 
with <j> m = 2-Km/Ngp (m = l,...,iV S p) as quadrature points to express 
the canonical trace in terms of grand-canonical traces. For example, the 
canonical trace of U a for A particles is given by 

1 Ns 

Tr A U a = — V e- l ^ A det (l + e^U,) . (7) 

y m— 1 

The multi-dimensional integral over the auxiliary fields in (J3J) is evalu- 
ated by Monte Carlo methods. The auxiliary fields are sampled according 
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to the distribution W a . Denoting the samples by {<Ji}, the expectation 
value in (j4|) is estimated from 



For a generic interaction, the sign $ CT can fluctuate from sample to 
sample and is in general a phase. At low temperatures, the fluctuations of 
the sign can become larger than its expectation value. This leads to large 
statistical errors in thermal observables and the breakdown of the method. 
When all components of the nuclear interaction discussed in Sec. 13.11 are 
attractive, we have Tr/To- > for any cr, and the interaction has a good 
Monte Carlo sign in the grand-canonical ensemble. Interactions with small 
bad-sign components can be treated by the method introduced in Ref. [22] 

3.3. Heat capacity 

In AFMC, we calculate the thermal energy E(T) = (H) as an observable 
and the heat capacity is determined from C = dE/dT. The statistical 
error in the numerical derivative can be reduced by an order of magnitude 
by calculating the energy at temperatures T and T + ST using the same 
Monte Carlo sampling and taking into account correlated errorsP^ To 
obtain the proper behavior of the heat capacity at higher temperatures, 
we use a method that combines correlated calculations in the truncated 
model space with independent-particle model calculations in the complete 
single-particle spaceP^ 

In Fig. Q] we show the heat capacities versus temperature for the even- 
even isotopes 56 Fe and 62 Fe (top left panel) and the odd-even isotope 57 Fe 
(bottom left panel). In the BCS approximation, the heat capacity is dis- 
continuous at the critical temperature (see the top right panel of Fig. [8] 
where we observe two discontinuities, one for the neutron pairing transi- 
tion and a second for the proton pairing transition). The AFMC results 
show significant suppression of the BCS heat capacity because of the large 
fluctuations in the pairing gap. However, in the even-even nuclei 56 Fe and 
62 Fe, there remains a shoulder in the heat capacity despite the large fluctu- 
ations. This shoulder structure is also refers to as an S-shape heat capacity 
and was observed experimentally in even-even rare-earth nuclei using the 
Oslo methocPS (see, e.g., the heat capacity of 172 Yb in the top right panel 
of Fig. [T]). For the odd-even nucleus 57 Fe (bottom left panel of Fig. [1}, the 
shoulder structure is suppressed, in qualitative agreement with the experi- 
mental result in 171 Yb (bottom right panel). 



(O) 
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Fig. 1. Heat capacity of even-even (top panels) and odd-even (bottom panels) nuclei. 
Top: the AFMC heat capacities of 56 Fe and 62 Fe (solid symbols) versus temperature T 
are compared with the heat capacity of 172 Yb extracted from experiments.EH Bottom: 
As for the top panels but for the odd-even nuclei 57 Fe and 171 Yb. The dashed lines are 
obtained in the independent particle model. The shoulder structure observed in the heat 
capacity of the even-even nuclei is a signature of pairing correlations. Notice that the 
shoulder structure is enhanced as neutron are added, turning it into a peak for 62p e [29] 



3.4. Spin distribution 

The spin distribution pj/p of nuclear energy levels {p and pj are, respec- 
tively, the total state density and density of levels with spin J) versus spin 
J at a given excitation energy E x can be calculated exactly in AFMC us- 
ing a spin projection method.^ We use the following identity for a scalar 
operator X 



TtjX = Tr M =jX - Tt m =,i+iX , 



(9) 
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where Trj denotes trace at a fixed spin J, while Tr^ denotes a trace at 
a fixed value M of the spin component J z . The projection Pm on a given 
value M of J z is accomplished using the Fourier sum 



Pm = 



1 



Jmax 

t E 

m— — J m 



e -i<j>mM e i<i>mJ z 



(10) 



where J max is the maximal value of the spin in the many-particle model 
space and 4> m = 27rm/(2J max + 1) are quadrature points. 

In Fig. [2] we show by solid squares the AFMC spin distributions for 
an odd-even nucleus ( 55 Fe), an even-even nucleus ( 56 Fe) and an odd-odd 
nucleus ( 60 Co) at excitation energies of E x — 4.39, 5.6, and 3.39 MeV, 
respectively.^ 
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Fig. 2. Spin distributions pj/p of energy levels in 55 Fe, 56 Fe and 60 Co. Solid squares 
are the AFMC results at the excitation energies indicated in the figure. The solid lines 
are empirical distributions (see text) deduced from global fits to experimentally known 
low-lying levels!^ The dashed lines are the empirical distributions but with the higher 
cr 2 values shown in the legends. These higher values are consistent with the higher 
excitation energies of the AFMC results. From Ref. 1331 



We compare our results with the spin-cutoff model obtained through 
the random coupling of the single- nucleon spins to total spin Jp3 In this 
model 

PJ (2J+1) _ZH+H . 
— = — t= — e , (11) 

p 2v / 2^cr 3 V ' 

where a is the spin-cutoff parameter. The AFMC spin distribution for 

odd-even and odd-odd nuclei are well described by the spin-cutoff model 
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(fTTj) with a fitted energy-dependent parameter a = u(E x ) (dashed lines in 
Fig. [2]). The energy dependence of a 2 (extracted from fits to the AFMC 
spin distributions) is shown in Fig. [31 The solid lines in Fig. [5] are the spin- 
cutoff formula with an empirical value of a 2 — 2.61A 028 as determined 
from global fits to spin distributions of experimentally known low-lying 
levels These distributions agree well with the AFMC distributions (solid 
squares) once the value of a 2 is scaled to larger values shown in the legends 
to take into account the higher excitation energies of the AFMC results. 
These empirical distributions with the scaled values of a 2 are shown by the 
dashed lines in Fig. [2J 
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Fig. 3. The squared spin-cutoff parameter <r 2 (top panels) extracted from the AFMC 
spin distributions (solid squares) and the J = pair correlations (At A) (bottom panels) 
versus excitation energy for 55 Fe, 56 Fe and 60 Co. The dashed lines correspond to cr 2 = 
IT IT? with rigid-body moment of inertia /. The open circles in the cr 2 panel of 55 Fe are 
experimental datafS] From Ref. I32I 



For even-even nuclei (e.g., 56 Fe), the spin-cutoff model works well only at 
higher excitation energies. As the excitation energy is lowered, an odd-even 
staggering effect as a function of spin is observed in the AFMC calculations. 
Such staggering effect was confirmed in the empirical studies of Ref. (33] and 
was parametrized by multiplying the spin-cutoff formula (|lll) by 1+a; where 
x = 0.227 (—0.277) for even (odd) non-zero spin values and x = 1.02 for 
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J = 0. This empirical distribution is shown by the solid line in Fig. [2] for 
56 Fe. The same distribution, scaled to a larger value of a 2 (dashed line) 
agrees well with the AFMC results. 

Fig. |3] shows the energy dependence of a as extracted from the AFMC 
spin distributions (solid squares). They are compared with the curves 
a 1 (Ex) calculated from a 2 — IT/ti 2 using a rigid-body moment of inertia / 
(dashed lines). For the odd-even and odd-nuclei we observe general agree- 
ment. However, for the even-even nucleus, we find that a 2 is suppressed at 
low excitation energies (below the pairing transition) when compared to its 
rigid-body value. The corresponding suppression of the moment of inertia 
is correlated with the onset of large neutron pair correlations (At A) (here 
A* = J2am a >o(~ i y a ~ maa ] a rn a a ] a -m a is the J — pair operator) at low 
excitation energies (see bottom middle panel of Fig. [3]) and is a signature 
of pairing correlations. 

3.5. Thermal moment of inertia 

The moment of inertia describes the response of the nucleus to rotations. At 
finite temperature and for a rotationally-invariant Hamiltonian, it is given 
by I = (3 (J 2 ). In Fig. [4] we show the AFMC results (circles) for / in iron 
nuclei. These results can be explained by a simple model (lines), in which 
we consider a monopole pairing Hamiltonian for a deformed nucleus.^ We 
then use the SPA together with a number-parity projection to describe the 
odd-even effects. In even-even nuclei / is suppressed at low temperatures 
because of pairing correlations, while in odd-even nuclei the suppression is 
weaker and in the limit T — ¥ we observe an enhancement because of the 
unpaired nucleon. 

4. Nano-scale metallic grains 

The spectra of nano-scale metallic grains were determined as a function 
of a Zeeman magnetic field by measuring the non-linear conductance of 
grains connected to external leads P In larger grains, a pairing gap was 
observed in the excitation spectrum of a grain with an even number of 
electrons. However, in smaller grains where the pairing gap A becomes 
comparable to the mean level spacing S, it is difficult to resolve such a gap. 
These grains describe the crossover between the bulk BCS limit and the 
fluctuation-dominated regime, in which BCS theory is no longer valid and 
the effects of pairing correlations are much more subtle. 
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Fig. 4. The thermal moment of inertia of iron isotopes versus temperature T. The 
symbols are the AFMC results and the lines are from a simple model described in the 
text. The dotted-dashed line are the rigid-body moments of inertia. Note the reentrant 
behavior of the moment of inertia of the odd-mass isotopes. From Ref. I36I 

4.1. Universal Hamiltonian 

Here we discuss a metallic grain whose single-electron dynamics is chaotic. 
The single-particle energies and orbital wave functions of such a grain ex- 
hibit mesoscopic fluctuations (from sample to sample or as a function of en- 
ergy) that follow random-matrix theory (RMT)f23 The two-body electron- 
electron interaction matrix elements, when expressed in the basis of single- 
particle eigenstates of the one-body Hamiltonian, fluctuate too. We can 
decompose these interaction matrix elements into an average and fluctuat- 
ing parts GMSS The fluctuating part of the interaction matrix elements can 
be shown to be suppressed by l/gr, where gr is the Thouless conductance 
of the grain. Here we consider the limit of large <?x, where the fluctuating 
part of the interaction can be ignored. 

The single-particle Hamiltonian together with the average part of the 
interaction defines the so-called universal Hamiltonian. 38 3 S This effective 
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Hamiltonian describes the low-energy physics of a chaotic metallic grains 
and contains three interaction terms: charging energy, spin exchange and 
pairing. For a fixed number of electrons, the charging energy is constant 
and the universal Hamiltonian has the following form 

H = E -gP^P- J S S 2 , (12) 

where 

p f = J24A ( 13 ) 

i 

is the pair operator and S is the total spin of the grain. 

Pairing correlations, which favors a superconducting minimal-spin 
ground state, compete with the exchange interaction, which favors a ferro- 
magnetic spin-polarized state. At zero temperature this competition leads 
to a coexistence regime in which the ground-state wave function is partly 
spin-polarized and partly pairedP^BI] Here we discuss signatures of this 
competition in thermodynamic observables (such as heat capacity and spin 
susceptibility) and their mesoscopic fluctuations. 

4.2. AFMC 

In the absence of exchange correlations (J s — 0), the universal Hamilto- 
nian (|12p reduces to a BCS-like Hamiltonian. Signatures of the BCS-like 
interaction in various thermodynamic observables of the grain were stud- 
ied by a number of methodsPSEUlSEl] Here we discuss an application of 
the AFMC method 4 ^ (an attractive BCS-like interaction has a good Monte 
Carlo sign in the density decomposition of the HS transformation). The 
band width of the grain is determined by the Debye frequency. However, in 
practical calculations we truncate the band width to 2N r + 1 single-particle 
levels, renormalizing the coupling constant g to keep the BCS gap fixed for 
the discrete system. For half filling of the band, the renormalized coupling 
constant g r is given b> * 46 * 47 ! 

9r 1 



arcsinh 



/ AV + l/2\ 

v ) 



(14) 



The heat capacity is calculated as in the nuclear case (see Sec. 13.31) . 
while the spin susceptibility \ = dM / (1B\b=o (M is a magnetization for a 
weak Zeeman field B) is calculated from 

X = 4/3/4 ((^ 2 ) - (S z ) 2 ) ■ (15) 
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Fig. 5. (a) Heat capacity C and (b) spin susceptibility x/XP (measured in units of the 
Pauli susceptibility \p = 2fi 2 B /S) of a metallic grain with A/5 = 3 and equally spaced 
single-particle spectrum. Solid (open) circles are the AFMC results for a grain with even 
(odd) particle number. The dotted lines are the results of the canonical independent- 
particle model and the dashed lines correspond to the BCS approximation. From Ref . 1461 



Here /j,b is the Bohr magneton and S z is the z component of the total spin. 

The AFMC heat capacity (left panel) and spin susceptibility (right 
panel) are shown in Fig. [5] for a grain with A/S = 3 and an equally spaced 
single-particle spectrum. The solid (open) circles correspond to an grain 
with and even (odd) number of electrons. The BCS heat capacity [dashed 
line in (a)] displays a discontinuity at the critical temperature for super- 
conductivity. In the finite grain, the heat capacity is a smooth function of 
temperature but exhibits an odd-even effect that is a signature of pairing 
correlations. In the odd grain, the spin susceptibility exhibits a re-entrance 
effect: it decreases as the temperature is lowered before diverging in the 
limit T -> 0. 



4.3. A finite-temperature method 

An attractive (ferromagnetic) exchange interaction leads to a sign problem 
in AFMC, but another quantum Monte Carlo method, suitable for a pairing 
interaction, is free of such a sign problem and was used to calculate thermo- 
dynamic observables of the grainP^ The energy eigenvalues of the universal 
Hamiltonian can also be determined by generalizing Richardson's solution^ 
to include the exchange interaction! 42 ! 50 ! However, to study the mesoscopic 
fluctuations, it is necessary to repeat the calculations for a large number of 
samples of the single-particle Hamiltonian. Both the quantum Monte Carlo 
method and Richardson's solution are computationally intensive. We used 
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instead a finite-temperature method,^ in which the exchange interaction 
is treated exactly by spin projection,^ and the pairing interaction is solved 
in the SPA+RPA together with a number-parity projection! 36 ' 521 ^ The 
number-parity projected SPA+RPA method was used in Ref. S3| to study 
odd-even effects in thermodynamic properties of a metallic grain in the 
absence of exchange correlations. 

4.3.1. Spin projection 

We treat the exchange interaction exactly using the spin projection method 
of Sec. 13.41 The partition function Z — Tr e~^ H of the universal Hamilto- 
nian at temperature T — 1//3 can be written as 

Z = J2 e^ s ( s+1 >Tr s e-^ BCS , (16) 
s 

where Trg is the trace over states with fixed spin S, and -Hbcs is the BCS- 
like pairing Hamiltonian 

H B cs= ]T eiclda-gPip. (17) 

i,er=t,4. 

Using the identity ([9]) with S replacing J and X = e _,9ifBCS , we can rewrite 
([TBI) in the form 



z = £ e ^s(s+i) {Zm=s Zm=s+1 ) . (18) 
s 

Here Zm = Tim ^e _/3i ^ BCS ^ with the trace taken at fixed value M of the 
spin component S z . The S z projection is given by a formula similar to (|10j) 
but with S z replacing J z . 

The spin susceptibility x can be calculated from Eq. (|15|) using (S z ) = 
(S 2 )/3 and (S z ) =0. We find 

X = ^Y, S ( S + 1 ^ 2S + i y JsSiS+1) ( Z M-S-Z M=s+1 ) . (19) 
6Zj s 

4.3.2. Number-parity projection 

We carry out particle-number projection in the saddle-point approxima- 
tion, where the canonical ensemble is approximated by the grand-canonical 
ensemble with an average particle number N. To describe the odd-even 
effects of pairing correlations, we use the number-parity projectioiPsEsHUl 

P v = \(l + e™*) , (20) 
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with r\ = 1 (r) = — 1) corresponding to projection on an even (odd) number 
of particles and N is the particle-number operator. 

4.3.3. Hubbard- Stratonovich transformation in the pairing channel 

The HS decomposition we use in the AFMC applications to nuclei and nano- 
scale metallic grains corresponds to a density decomposition, in which the 
auxiliary fields cr(r) are densities. The reason for doing so is that there is no 
Monte Carlo sign problem in such a decomposition. Here we use a pairing 
decomposition, in which the auxiliary field is the complex pairing field A(r). 
In this pairing decomposition, the propagator of the BCS Hamiltonian in 
the grand-canonical formalism is given by 

f> 

n/,r A\ P - fdr\A(T)\ 2 /g 

e -0(H BCS -»N) = / P [ A;A *] e o /9 Ua , (21) 

- / dr -Ha(t) 

Here J7a = Te (T denotes time ordering) is the propagator for 

the one-body Hamiltonian 

#a(t) = [( e * - ^ - f ) ( C *W + C W) 

i 

- A(r) C t T ^-A*(r)c 4 c it + |] . (22) 

4.3.4. SPA + RPA 

The number-parity and ^-projected partition function Z 7h M — 
Tr 



can be calculated using the HS transformation 
([2~Tj) as follows. We expand the pairing field in a Fourier series A(r) = 
Ao + ^ r ^ A r e" JrT , where uj r = 2jrr/f3 (r integer) are bosonic Matsubara 
frequencies. For each static fluctuation Ao, we perform the integration over 
A r (r 0) in the saddle-point approximation and then keep the exact inte- 
gral over |Ao| (the integration over the phase of Ao is trivially done). The 
projection on particle number N is performed in a saddle-point approx- 
imation, leading to the following expression for the ./V-particle partition 
function at fixed r\ and 



Z 



N,n,M 





/2tt 


d 2 F 




1 /5d|A | 2 






) 



-1/2 



x e -(^/ 9 )|Aol 2 e -0^ Zr)MAo) c^(A ) . (23) 
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Here Z v m(Ao) i s the number-parity and ^-projected partition function 
for a static fluctuation Ao 

-i<j> m M 




-fi—Ei) 



E 



2(25 

max 



i) 



(24) 



where Et 



ti — \i — |) +|A | 2 are the quasiparticle energies for the 
given Ao- The factor C?j^(Ao) is a local RPA correction factor, arising 
from the small-amplitude time-dependent fluctuations 



(Ao) 



n 



sinh(/3£; j ) 



2Ei 



sinh 



2 



(25) 



The local RPA frequencies ±fij are the eigenvalues of the 2i\T sp x 27V sp RPA 
matrix (N sp is the number of single-particle orbitals) 



2Ei5ij 



|/i(7i7j 
1) 



with 7i = (a - fj, 



jfi {lilj 

f ) /Ei and 



\fihi7j + 1) 



2EiSij 



fi = 



1 ainZ„, M (Ao) 



P BE, 

F in Eq. (|23|) is the grand-canonical free energy 
|Ao| 2 



G 



^2 (e i -M)-2^ 1 ln^2cosh^^ 



(26) 



(27) 



(28) 



and /x is the chemical potential determined (for each Ao) from the particle 
number equation N = —dF/dfi. 

The partition function and spin susceptibility for a metallic grain with 
N electrons and number parity r\ are calculated from Eqs. (TT81 and (|19p . 
respectively, where Zm is the partition function Zjv.tj.m m Eq. (1231) . 

The SPA+RPA method breaks down below a critical temperature when 
a fluctuation A r for a certain value of Ao becomes unstable. A method was 
recently proposed^ to overcome this problem. 

In Fig. [5] we compare the number-parity projected SPA+RPA re- 
sults (symbols) with the exact results obtained from Richardson's solution 
(lines). We see that this approximation is very good except that it breaks 
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down below a certain critical temperature. However, the odd-even signa- 
tures of interest are observed above this temperature. In using Richardson's 
solution we calculate all the eigenvalues of the universal Hamiltonian below 
an energy cutoff of ~ 30 8 so the corresponding results in Fig. |6] are no 
longer accurate above T/8 ~ 1.5. The correct heat capacity is then given 
by the number-parity projected SPA+RPA results. 




Fig. 6. The results of the spin and number-parity projected SPA+RPA method (solid 
circles for an even particle number and solid squares for an odd particle number) are 
compared with exact results obtained from Richardson's solution (solid lines for an even 
grain and dashed lines for an odd grain): (a) heat capacity for a grain with A/5 = 3 
and J s /S = 0.5; (b) spin susceptibility x/XP f° r a grain with A/5 = 0.5 and J s /5 = 0.5. 
The single-particle spectra correspond to specific RMT realizations. From Ref. 1511 



4.4. Thermal observables: heat capacity and spin suscepti- 
bility 

We used the method of Sec. l4.3l to study the mesoscopic fluctuations of the 
heat capacity and spin susceptibility of a metallic grain for a large num- 
ber of realizations of the single-particle RMT spectrum.^ The results are 
summarized in Fig. [7] for both even and odd grains and for different values 
of A/ 8 and J s /8. The symbols are the average values over the ensemble 
and the vertical bars describe the standard deviations of the corresponding 
quantities. The lines are the results for an equally spaced single-particle 
spectrum. 

As already discussed in Sec. 14.21 pairing correlations lead to odd-even 
effects in the heat capacity and spin susceptibility. Here we study how ex- 
change correlations and mesoscopic fluctuations affect these number-parity 
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A/5 =1.0 A/5 = 3.0 A/5 =1.0 A/5 = 3.0 




T/5 T/5 

Fig. 7. Mesoscopic fluctuations of (a) heat capacity C and (b) spin susceptibility x/Xp 
in metallic grains with A/S = 1 (left columns) and A/S = 3 (right columns). Results 
are shown for different value of J s /S and for both even grains (circles) and odd grains 
(squares). The symbols and vertical bars describe, respectively, the average values and 
standard deviations (over the ensemble of single- particle spectra). The lines are the 
results for an equally spaced single-particle spectrum and the dotted-dashed lines are 
the grand-canonical BCS results. From Ref. 1511 

dependent signatures of pairing correlations. 

In general, the exchange interaction shifts the odd-even effects in 
the heat capacity and spin susceptibility to lower temperatures. In the 
fluctuation-dominated regime A/S < 1, exchange correlations suppress the 
odd-even effect in the heat capacity as well as the reentrant behavior of the 
spin susceptibility for an odd number of electrons. However, for A > S, ex- 
change correlations enhance the shoulder in the even particle-number heat 
capacity and can turn it into a peak (see the panel in Fig. [7] with A/S = 3 
and J s /S = 0.8). Similarly for A > 5, exchange correlations enhance the 
reentrant effect for an odd number of electrons (see the panel with A/S = 3 
and J s /S = 0.6). 

In the fluctuation-dominated regime, A/ 5 < I, the mesoscopic fluctu- 
ations of the heat capacity can wash out the odd-even effect for moderate 
strengths of the exchange interaction. In this regime we also observe large 
fluctuations of the spin susceptibility as the exchange coupling constant 
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increases. 



5. Conclusion 




Fig. 8. Thermal signatures of pairing correlations in nano-scale metallic grains with 
equally spaced single-particle spectrum and A/5 = 1 (left panels), and in iron nuclei 
(right panels). Top panels: the heat capacity. Bottom panels: the spin susceptibility 
(bottom left for a metallic grain) and the thermal moment of inertia (bottom right for 
iron nuclei). The symbols are AFMC calculations. The solid lines in the right panels 
(metallic grains) are calculated from Richardson's solution. The solid lines in the bottom 
right panel (nuclei) are from the simple model of Ref. 36 The dashed lines are the results 
of the BCS approximation. 

We discussed thermal signatures of pairing correlation in nuclei and 
nano-scale metallic grains in the crossover between the bulk BCS limit and 
the fluctuation-dominated regime where BCS theory breaks down. We used 
methods that go beyond the mean-field approximation: AFMC and a finite- 
temperature method that takes into account thermal and small-amplitude 
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quantal fluctuations together with spin and number-parity projections. A 
summary of some of our results is shown in Fig. [8l in which we compare 
thermodynamic properties of a metallic grain with A/8 — 1 to similar ther- 
mal properties in iron nuclei. An odd-even effect in particle number is 
observed in the heat capacity of both the metallic grain and the nucleus. 
Similarly, we compare the spin susceptibility (response to an external mag- 
netic Zeeman field) of metallic grains with the moment of inertia (response 
to rotations) of nuclei. In both metallic grains and nuclei, we observe a 
reentrant effect for an odd number of particles. 

The results in Fig. [5] also demonstrate the large discrepancy between 
the BCS approximation and the exact results for various thermodynamic 
observables, emphasizing the necessity to use methods beyond the mean- 
field approximation. 
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